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Abstract. We present a study of the medium distant cluster of galaxies Abell 1835 based on Xmm-newton 
data. The high quality of Xmm-newton data enable us to perform spectro-imaging of the cluster up to large 
radii. We determine the gas and total mass profiles based on the hydrostatic approach using the /3-model and the 
temperature profile. For the determination of the temperature profile of the ICM, which is needed for the mass 
determination, we apply a double background subtraction, which accounts for the various kinds of background 
present (particle and astrophysical background). We find a basically flat temperature profile up to 0.75 r2oo with a 
temperature decrease towards the center linked to the cooling flow. We obtain a gas mass fraction of (20.7±3.7)%, 
which is a lower limit on the baryon fraction in this cluster. Using this value as baryon fraction for the entire 
universe, we obtain by combining our results with results based on primordial nucleosynthesis, an upper limit for 

— 1/2 

Sim < 0-5h^Q , which is in good agreement with other recent studies. 
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' 1. Introduction 

. Clusters of galaxies are dark matter dominated and thus 
' are ideal objects to obtain information on this kind of 
material. A straightforward way to determine the mass 
profile of a cluster, which consists of roughly 80% dark 
matter is to use the hydrostatic approach based on the 
hot x-ray emitting intra-cluster medium (hereafter icm). 
In order to get reliable results with this approach it is 
necessary that the examined galaxy cluster is in a relaxed 
state without ongoing merger activity, which might cre- 
ate shock waves and bias the mass estimate. The validity 
of the hydrostatic approach for relaxed clusters has been 
tested via hydrodynamic simulations (Schindler 2001 ; 
Evrard, Metzler & Navarro 1996). To solve the hydro- 
static equation one needs to know the density as well as 
the temperature distribution of the ICM. The latter one 
requires x-ray observatories, which provide at the same 
time good spatial as well as spectral resolution. X-ray ob- 
servatories such as ROSAT, ASCA, and Beppo-Sax only 
provided limited combinations of these two requirements 
and thus the determination of cluster temperature profiles 
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based on these data were subject to large error bars and 
sometimes different studies on the same object with the 
same data came to very different conclusions (Irwin & 
Bregman 2000 ; Markevitch et al. 1998). A subject of hot 
debate was the question whether there exists a univer- 
sal declining temperature profile in clusters as suggested 
by Markevitch et al. (1998) or not. Recently launched x- 
ray satellites namely Chandra (Weisskopf et al. 2000) and 
Xmm-newton (Jansen et al. 2001) fuUfill now aU require- 
ments for an accurate determination of the temperature 
distribution in clusters up to large radii. 

We present here a spectro-imaging study and subse- 
quent mass determination analysis of the medium distant 
(z=0.25) relaxed galaxy cluster Abell 1835 (Allen et al. 
1992) observed during the Xmm-newton Performance 
Verification phase. Abell 1835 is a relaxed cluster, which 
hosts a cooling flow in its center (Allen et al. 1996). It 
belongs to the cooling flow clusters for which no multi- 
phase gas down to low temperatures has been found in 
the Xmm-newton rgs data (Peterson et al. 2001), and 
which casts serious doubts on current cooling flow models 
(Fabian et al. 2001, Molendi & Pizzolato 2001). 

We determine the total and gas mass profile of the 
cluster up to 1.7Mpc which corresponds to 0.75 virial ra- 
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dius (see section ^ . The ratio of gas mass over total mass 
allows us to determine a lower limit for the baryon fraction 
of the cluster which can be considered as representative of 
the mean baryon fraction of the Universe (see White et al. 
1993 and references therein). Knowing the baryon density 
of the Universe based on studies of primordial nucleosyn- 
thesis we can give constraints on the cosmological matter 
density parameter fl^. 

Our paper is organized as follows : in section |2| we de- 
scribe the observation. In section || wc present the data 
analysis with specific emphasis on vignetting corrections 
and flare rejection. In section ^, we present first results and 
a method to correct for different background components : 
a) particle induced background and b) astrophysical back- 
ground, which varies across the sky (Snowden et al. 1997). 
The method is based on a double background subtraction. 
Section ^ presents our gas and total mass determination 
and estimate of the baryon fraction of Abell 1835. This is 
followed in section |6| by the discussion of our results and 
our conclusions in section 0. 

Throughout this article, we use a cosmology with Ho = 
50km/s/Mpc, = 1 and fl\ = (qo — 0.5). One arc- 
minute corresponds to 297 kpc at the cluster redshift of 
0.25. 

2. Observation 

Xmm-newton observed AbeU 1835 in revolution 101 on 
June 28th 2000. The observation (id 98010101) was taken 
during the phase of Performance Verification. The total 
exposure time was 60ksec. Since we are interested in spa- 
tially resolved spectroscopy we concentrate here on data 
from the EPiC-instruments (Turner et al. 2001 ; Striider 
et al. 2001). During the observation of Abell 1835 the 
MOSl camera was operating in Large Window mode, which 
shows a relatively high internal background level, and for 
which we cannot yet properly account for. Since our work 
is very sensitive to background uncertainties we discard 
MOSl data of Abell 1835 in the following and only use 
data coming from the MOS2 and pn camera. For all cam- 
eras thinI filter was used. 

3. Data processing 

We generate calibrated event files using the tasks emchain 
and epchain in the XMM SAS version 5.0 0. For the MOS2 
data set, we take into account event patterns to 12 and 
for pn data only single events (pattern ==0). 

3.1. Flare rejection 

The light curve of this observation (see figure |^) is not con- 
stant and large variations of intensity are visible. These 
variations, which we will call "flares" in the following, are 
caused by soft energy protons produced by solar activity, 

^ http://xmm.vilspa.esa.es/sas 
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Fig. 1. Light curve derived from MOS2 camera and event 
energy between 0.3 and 12 keV. 
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Fig. 2. Effective exposure time versus count rate thresh- 
old for several MOS data sets. For the Lockman Hole ob- 
servations, the revolution 70, respectively 71 was used for 
MOSl, respectively MOS2. The count rate threshold was 
defined as the number of events per 100 s time intervals in 
the 10 to 12 keV energy band. 



and which are only visible outside the Earth's magneto- 
sphere, which corresponds to the orbit of Xmm-newton 
(see Ehle et al. 2001). 

Unfortunately, the background induced by flares can- 
not be easily accounted for since it shows spectral vari- 
ability. Furthermore, the flares reduce the observed signal- 
to-noise ratio. Therefore it is best to discard flare periods. 
This reduces the effective exposure time, but improves the 
signal-to-noise ratio significantly. 

For our analysis we only consider events inside the 
field-of-view (fov) and ignore all events outside the FOV. 
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Fig. 3. Same as figure ^ except it is for pn data sets and 
that the threshold was defined as the number of events 
per 100 s time intervals in the 12 to 14keV energy band. 



Flares show a hard component observable at high en- 
ergies easily detectable above lOkeV where the mos cam- 
eras are hardly sensitive to X-ray photons. We choose the 
10 to 12keV energy band to monitor the background for 
the MOS2 camera. In order to optimize the background 
rejection, we analyze the relationship between effective 
exposure time and the threshold we apply. The results 
are shown in figure ^ for Abell 1835 and Lockman Hole 
(revolution 70 and 71) observations. We can see that the 
behaviour of exposure time versus rejection threshold is 
similar for all observations, with a flattening at a thresh- 
old of about 15 counts/lOOs. Thus we can say that this is 
the optimal threshold for flare rejection, since it limits ex- 
posure time loss. A higher threshold does not significantly 
increase the exposure time and a lower threshold limits 
severely the remaining effective exposure time. Therefore 
we adopt in the following analysis a threshold of 15 
counts/lOOs in the 10 to 12keV energy band. The same 
procedure is applied for the pn data (see figure^). The pn 
camera provides higher quantum efficiency at high ener- 
gies than the MOS detectors. Because of this we shift our 
energy window for monitoring the background level to the 
12 to 14keV energy band. Our adopted threshold is 22 
counts/ 100s. 

Of course, the adopted threshold criterion for flare re- 
jection must also be applied to the observations chosen to 
determine the background. 

After applying this flare rejection method, the remain- 
ing exposure times are 25750 s for MOS2 and 25950 s for 
pn. 




Fig. 4. Mos2-f pn vignetting corrected image of Abell 
1835 in sky coordinates in the 0.3 to 3keV energy band. 
This image is filtered with a cr=2.2' Gaussian. The con- 
tours come from an 0.3 to 3keV only MOS2 image where 
the gaps between CCDs cannot change the shape of these 
contours. The contour levels are 0.25, 1, 2, 4, 8, 16 and 28 
photons by 1.1" xl.l" pixel. 



3.2. Vignetting correction 

The collecting area of the x-ray telescopes aboard Xmm- 
newton is not constant across the field of view. The ef- 
fective area decreases with increasing off-axis angle and 
is furthermore photon energy dependent. This effect (vi- 
gnetting) must be taken into account in the case of ex- 
tended sources, like galaxy clusters, which cover large frac- 
tions of the field of view. 

To correct for this vignetting effect, we use the method 
described by Arnaud et al. (2001a). Briefly, this method 
consists of calculating a weight parameter for each event. 
According to event energy, this parameter is defined as the 
ratio of on-axis effective area to effective area at the event 
position on the detector. Thus, all weight parameters are 
greater or equal to unity. 



4. Data analysis 

4.1. IVIorphology 

As can be seen in figure ^, Abell 1835 is a compact clus- 
ter. There is no strong evidence of substructures which 
suggests that the cluster is in a relaxed state. The inten- 
sity is very peaked in the center due to the presence of a 
cooling flow (Allen et al. 1996; Schmidt, Allen & Fabian 
2001; Peterson et al. 2001). 
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4.2. Spectra 

Since we apply the described vignetting correc- 
tion method, we can use the on-axis response 
(i.e. effective area x redistribution matrix) files : 
m2_thinlv9q20t5r6_alL15.rsp for the MOS2 camera and 
epn_fF20_sY9_thin.rsp for the pn. There exist different re- 
sponse matrix files for the pn-chips, which correspond to 
the different nodes. We use Y9, which is the node farthest 
away from the readout, with lowest resolution. We tested 
the different response files of the different nodes and 
could not find any detectable difference when performing 
spectral fits. Due to the remaining uncertainties of the 
response files at very low energy, we exclude events with 
energy below 0.3 keV for all spectral fits. 

4.2.1. Background treatment 

Due to the fact that clusters of galaxies are extended and 
generally low surface brightness sources, the correct deter- 
mination of the background is important for data fitting 
and modelling. Since cluster emission decreases with ra- 
dius from the center, the background component becomes 
more and more important with increasing radius. 

After removing the flare component by excising high 
background intensity time intervals, two other background 
components remain : 

— high energy particles like cosmic-rays pass through 
the satellite and deposit a fraction of their energy 
on the detector. They are not affected by the tele- 
scope vignetting. These particles induce instrumen- 
tal fluorescent lines coming from material aboard the 
satellite. This kind of background is well monitored by 
blank-sky observations based on several high galactic 
latitude observations in which x-ray sources were ex- 
cised. These observations were compiled by D. Lumb 
(2002)1 for each epic camera ; 

— the astrophysical x-ray background. The soft com- 
ponent of this background is position dependent and 
changes across the sky. Therefore it is not necessar- 
ily well monitored by blank-field observations taken at 
different positions on the sky. 



A description of the following method to correct 
for each background component can also be found in 
Majerowicz & Neumann (2001) and Pratt et al. (2001). 
A similar method was used by Markevitch & Vikhlinin 
(2001) for Chandra data. 

In the data sets of Abell 1835 and of the blank-sky 
observations, which we use for background subtraction, 
we attribute to each event a weighting factor to correct 
for vignetting (see section 3.2). This is obviously wrong 
in case of particle events, which are not vignetted. If the 
particle background (after flare rejection, see section 3.1 ) 
is constant with time in terms of spectrum shape, what 
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Fig. 5. A MOS2 (lower) and pn (upper) combined fit of 
fully background corrected spectra in the outside region 
of the cooling flow (see section 4.4). 



^ They can be retrieved from the Vilspa ftp site : 
ftp://xmm.vilspa.esa.es/pub/ccf/constituents/extras/background 



we assume here, we can correct for the falsely attributed 
weighting factor by subtracting the background spectrum 
in the same detector region as the source spectrum. In or- 
der to account for possible long time intensity variations in 
the particle induced background, we determine the count 
rates of source and background observation at high ener- 
gies (> 10 keV), where the quantum efficiency for X-ray 
photons are low. The ratio of the countrates is used to 
normalize the corresponding background in the spectral 
fits. For the MOS2 we use the 10-12 keV band for deter- 
mining the background intensity, and for the pn-data the 
12-14 keV band, respectively. In our case, a normalization 
factor of 1 is found. With this method we correct for the 
particle background. However we do not yet necessarily 
account for the spatial variations of the CXB. 

Abell 1835 is located in a peculiar region of the sky, 
where an excess of low energy CXB can be seen (Snowden 
et al. 1997). To correct for this local excess of the CXB, we 
use regions of the detector in which x-ray emission from 
the cluster is negligible. We assume in the following that 
the CXB does not change significantly across the 30' of 
detector field of view. We extract a spectrum in the annu- 
lus between 7' and 12' from the cluster center. From this 
spectrum outside of the cluster, we subtract the blank- 
sky spectrum extracted from the same detector region. 
The residual spectrum represents the difference of the 
CXB between the observation of Abell 1835 and the blank- 
sky observations. These residuals were subsequently sub- 
tracted from the cluster spectrum, from which we already 
subtracted the blank-sky spectrum in the same detector 
regions. We thus yield spectra which are now properly 
corrected for the two background components described 
above. In figure H, we present background corrected spec- 



S. Majerowicz et al.: Xmm-newton observation of Abell 1835. 



5 



tra from MOS2 and pn data wliicli will be used in section 



4.4 



for the estimate of the mean temperature outside the 
cooling flow region. 

For the pn data there exists another source of con- 
tamination, the so called out of time (oot) events, i.e. 
events counted during the read-out (see Striider et al. 
2001). They are visible as a jet -like structure in figure ^. 
To estimate the impact of these OOT events on the temper- 
ature determination, we extract spectra in regions defined 
in table |l] in which we exclude the OOT affected regions. 
Comparing the fit results of those OOT corrected spec- 
tra with the non-corrected OOT spectra we find a slight 
increase in the temperature estimates of a few percent, 
however, always well within the error bars. We thus con- 
clude that the OOt's do not play an important role in the 
temperature determination for this particular observation. 
In the following analysis of pn data, we use the non-OOT 
corrected spectra. 



1 ooc 

Radius (kpc) 



4.2.2. Temperature profile 

Abell 1835 appears to be a relaxed cluster of galaxies (see 
section 4.1). Therefore we assume that the temperature 



structure of this cluster is spherically symmetric. We ex- 
tract cluster spectra in concentric annuli around the clus- 
ter center. The size we choose for the annuli is a compro- 
mise between number of source counts to have sufficient 
statistics on one hand and small enough regions to deter- 
mine the temperature profile with small bins on the other 
hand. In our analysis, we exclude the most luminous point 
sources which are located inside the field of view of the 
cameras. 

We group spectra such that each bin has a signal-to- 
noise ratio greater than 3cr after background subtraction. 
This grouping allows us to assume Gauss statistics, which 
is important for fitting. 

For the fit to the spectra, we use xspec version 11.0.1 
(Arnaud 1996) and we model the obtained spectra with : 



Model = Wabs(nH) x Mekal(T, A) 



(1) 



where Wabs is a photo electric absorption model 
(Morrisson & McCammon 1983) and Mekal a single tem- 
perature plasma emission model (Mewe et al. 1985 ; 
Kaastra 1992 ; Liedahl et al. 1995). The free parameters 
are the galactic hydrogen column density value uh, the 
x-ray temperature kT and the metallicity (abundance) 
A. 

First we fit all the spectra with a fixed galactic absorp- 
tion uh = 2.24x 10^° cm~^ given by Dickey & Lockman 
(1990) from 21cm line measurements. Results are shown 
in table |l]. In a second spectral analysis we leave the nn 
value as a free fit parameter (see table |^). The resulting 
fit temperatures agree very well with those obtained with 
fixed uh which gives us high confidence in our spectral 
analysis and our background subtraction. Figure ^ shows 
the fit results for uh which are in good agreement with 
the measurements from the 21 cm line. Figure ^ shows the 
resulting temperature profile up to a physical radius of 1.8 



Fig. 6. Abell 1835 temperature profile (see table |i|) de- 
rived from a combined fit of the MOS2 and pn cameras 
where the nn value is equal to the galactic value at the 
sky position of Abell 1835. The mean temperature (see 
section 4.4) is represented by the dashed line. 
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Fig. 7. Abundance profile (see table |l]) and its mean value 
outside the cooling flow (dashed line, see section 4.4) in 
units of solar metallicity up to 1 Mpc. 

Mpc derived from the table 0. Figure |^ shows the metal 
abundance profile up to 1 Mpc. 

4.3. PSF corrected temperature profile 

Abell 1835 is a medium distant galaxy cluster with a very 
peaked surface brightness profile (see figure |l^) due to the 
presence of a cooling flow in the centre. Because of this, the 
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Annulus (') 

ri r-z 



MOS2 



T (keV) ixi^Jdof) 
pn 



MOS2+pn 



A (solar) 
MOS2+pn 



0.0 
0.25 
0.75 

1.5 
2.25 
3.33 



0.25 
0.75 
1.5 
2.25 
3.33 
6.0 



4.81 ±0.25 (1.17/307) 
6.67 ±0.29 (0.98/371) 
7.29lo 39 (0.95/325) 
7.24+E!j3 (0.91/229) 
6.87; 
6.98 



^■^^ (1.10/164) 
,.S9 (0.95/53) 



-0.95 
+3.35 



4.43 ± 0.12 (1.17/706) 
6.16 ±0.20 (1.12/933) 
6.76 ±0.27 (1.02/749) 
7.15l»i° (1.02/509) 
8.3lj;i:J^ (0.90/339) 
7.80tigg (1.15/84) 



4.52 ±0.10 (1.18/1014) 
6.31 ±0.16 (1.09/1305) 
6.94 ±0.22 (1.00/1075) 
7.19 ±0.38 (0.99/739) 
7.68t^;?^ (0.97/504) 
7-3lt?:33 (1.06/138) 



0.42 ± 0.05 
0.30 ± 0.04 
0.23 ± 0.05 
0.32 ± 0.08 
0.27 ±0.13 
< 0.9 



Table 1. Results for the isothermal model fits derived from equ. (0) where the free parameters are the temperature T 
and the abundance A. The galactic hydrogen column density value nn is frozen to the galactic value of 2.24x 10^*^ cm~^ 
(Dickey & Lockman 1990). The spectra are extracted in annuli with inner radius ri and outer radius r2. The errors 
have a confidence level of 90%. 
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Fig. 8. Absorption column density nn profile from ta- 
ble ||. The dashed line represents the galactic value of 
2.24xl020cm-2 (Dickey & Lockman 1990). 



Annulus 


T (keV) 


nn (10^" cm") 




1 


^•"'-'-0.12 


1.90 ±0.30 


1.17/1013 


2 


6.61 ±0.20 


1.64 ±0.23 


1.07/1304 


3 


r, o, +0.33 
' •■J^-0.29 


1.56 ±0.32 


0.99/1074 


4 
5 
6 


V Q/1+0.66 
'•■J*-0.47 

7 c;8+i-i3 

'•JO_0.96 


1.95 ±0.52 

^•J^"_0.78 
o qn+2.3S 
^•'^"-l.SO 


0.99/738 
0.97/503 
1.07/137 



Table 2. Spectral fit results with nn value as free fit pa- 
rameter (MOS2±pn combined fits). 



point spread function (psf) of Xmm-NEWTON (Griffiths 
et al. 2002, Ghizzardi 2001) potentially blurs the central 
regions of the cluster, which might influence the observed 
temperature profile (see also Markevitch 2002) . In order to 
assess this influence we perform another spectral analysis, 
in which we take into account the PSF of Xmm-newton. 
We apply the method which is extensively described by 
Pratt & Arnaud (2002). 



For our analysis we assume that the PSF is constant 
across the extent of the cluster and that its shape is en- 
ergy independent. We calculate the flux redistribution of 
all temperature bins due to blurring. This computation is 
based on a convolution of a double /3-model, which is fitted 
to the surface brightness profile and which gave the follow- 
ing parameters (this fit was performed taking the PSF into 
account). Subsequently, based on these model parameters 
(An = 0.49, r,,i„ = O.r, /3out = 0.74, r^out = 0.79') we cal- 
culate the flux contribution of bin i to bin j. The result 
can thus be interpreted as a matrix. The corresponding 
redistribution elements are shown in figure ^ for the pn 
data. The redistribution for MOS2 is similar. As an illus- 
tration of the importance of the PSF effect in the centre : 
only 65% of the flux of the second bin (i = 2) are actually 
observed in the second bin (j — 2). The remaining 35% 
are distributed in other bins [i = 2,j ^ 2). 

After the determination of the redistibution or matrix 
elements wc fit simultaneously the different spectra of our 
six temperature bins. We take into account the different 
flux contributions of each bin by fixing the ratios of nor- 
malization of the fitted spectral models based on our cal- 
culated redistribution matrix. Our resulting Xred ^o^' 
combined fit is ~x^^^ = 1.05 for 4775 degrees of freedom. 
The resulting temperature profile is displayed in figure ^ 
as well as in table ^. Briefly, when we take into account 
the PSF we obtain a steeper temperature gradient in the 
cooling flow region of the cluster. The temperature esti- 
mates at larger radii are hardly affected by the PSF and 
are within the error bars of the non-PSF corrected analy- 
sis. Our study presented here is focussed on the external 
regions of the cluster, therefore we conclude that the PSF 
does not have an important impact on our analysis (see 
below). 

In the following analysis, we will exclusively use the 
PSF corrected temperature profile. 

4.4. Mean temperature outside the cooling flow 

We use the background subtracted spectra for MOS2 and 
pn data in an annulus between 1' and 6' — between 
300 and 1800 kpc — around the cluster center to obtain 
the mean temperature of the cluster outside the cool- 
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Annulus T (keV) Abundance (solar) 



1 


4.11 ± 0.12 


0.43 ± 0.05 


2 


8.03 ± 0.39 


r, 00+0.06 
""J^-O.OS 


3 


7.12 ±0.36 


0.19 ±0.08 


4 


'•'^-0.77 


0.37t°:l^ 


5 


o cc+i.ae 

o-'JJ_1.84 


< 0.5 


6 




< 1.0 



Table 3. Temperature estimate in each annulus for PSF 
corrected MOS2+pn combined fit. 



0.0 0.5 1.0 l.b 2.0 2,5 3,0 

Radius (kpc) 

Fig. 9. Redistribution of the cluster flux for the pn camera 
in the radial bins we use to derive the temperature profile. 
The triangles represents the flux coming from the bin, the 
squares from all inner bins to the bin, the losanges from 
all outer bins to the bin and the crosses from the central 
bin to the bin. 



:coo 

Radius (kpc) 



Fig. 10. Temperature profile of Abell 1835 corrected for 
the PSF (from table |[ squares) comparing to the previous 
one (from table |^, crosses). The dashed line represents 
the mean temperature outside the cooling flow region (see 
section I 



ing flow region. We perform a combined fit using the 
model described above in equ. (jl|) with a uh value fixed to 
2.24xl0^°cm"2. We find a temperature of 7.6 ± 0.4keV 
and an abundance of 0.25±0.05 with 90% confidence level. 



4.5. Surface brightness profile 

We extract surface brightness profiles from the cluster 
data. In order to account for vignetting effects we use 
again the photon weighting method described in section 
|3.2| . In our specific case using the weight factor for each 
event to account for vignetting, the count rate CRi in the 
i-th annulus is defined as : 



CR, 



1 



^cxpo -^i 



(2) 



where n is the number of events in the i-th annulus, Wj 
the event weight parameter, Ai the surface of the selected 
annulus and Toxpo the exposure time. In this case, the 
Poisson error becomes : 



expo 



A, 



(3) 



Serendipitous sources in the field of view are excluded. 
Figure |ll| shows the surface brightness profile of Abell 
1835 derived from the MOS2 and pn cameras in the 0.3 to 
3 keV energy band. 

The surface brightness profile is fitted with a /3-model 
(Cavaliere & Fusco-Femiano 1976) in which the surface 
brightness S is defined as : 




-3/34 



B 



(4) 



where 5*0 is the central intensity, the core radius, [3 the 
slope parameter and B the background intensity. All these 
parameters must be fitted. 

Abell 1835 hosts a cooling flow in its center (see above). 
Figure ^ shows the characteristic peak of intensity for 
cooling flows towards the center (e.g. Fabian 1994). In or- 
der to avoid fitting in the cooling flow region we choose 
to exclude a region with a radius of roughly 0.1r2oo (see 
equ. (||)) — i.e. 0.7' — around the center for the /3-model 
fit. Indeed Neumann & Arnaud (1999) showed that be- 
low 0.1 r2oo the surface brightness profiles for nearby clus- 
ters show a large dispersion due to cooling flow but they 
look remarkably similar above 0.1r2oo- The best fit pa- 
rameters are (3 = 0.704 ± 0.005 and a core radius — 
202.3 ± 7.1 kpc (i.e. = 0.681 ± 0.024') for the MOS2- 
data. The uncertainties are 68% confidence level. We also 
perform a /3-model fit with the surface brightness profile 
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Fig. 11. Abell 1835 surface brightness profiles from MOS2 
(upper panel) and from pn (lower panel) in 0.3 to 3keV 
energy band (I'corresponds to 297 kpc). The error bars are 
one a based on Poisson statistics. The resulting /3-model 
fits are shown in full line. 
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Fig. 12. Abell 1835 gas mass profile. The error bars (dot- 
ted line) are a confidence level of 68%. The gas mass profile 
(full line) is derived from the best fit parameters. 



extracted from the pn data for consistency. We find in 
this second case /? = O.TOSt^ J^Jy, = 196.9^^:3 kpc (i-^- 
Tc = 0.663lo^Qig '). These fit parameters are thus in good 
agreement with those found with the MOS2 camera. 

Our found parameters of /3 and rc are in relatively good 
agreement with the results found by Neumann & Arnaud 
(1999) for nearby clusters. 



5. Mass analysis 

We assume in the following that the cluster is virialized up 
to r2oo , the radius in which the mean density of the cluster 
is 200 times the critical mass density. For the calculation 
of ^200 — ^viriai wc usc the relation by Evrard et al. (1996) : 



7-200 = 3690 



T 



lOkeV 



(1 + ^)" 



(5) 



Using the mean temperature of 7.6 keV (see section , 
we find r2oo = 2293 kpc. 

5.1. Intra cluster medium mass 

Assuming that the radial ICM distribution follows a f3- 
model, the electron density rig can be written : 



(6) 



From the electron density we calculate the total gas den- 
sity of the ICM and determine the ICM mass by inte- 
grating equ. (H). We obtain a central electron density 




Ueo = i.4(' ± u.uz X iu "cm " Confidence level) as- 

suming a temperature of 7.6 keV and using the /?-model 
parameters determined in section 4.5. Temperature varia- 



tions in the uncertainties do not change our values for the 
electron density. The error bars of n^Q are calculated from 
the errors of the /3-model parameters. We verify that the 
cluster temperature has no effect on the estimate of Ueo 
with the selected energy range. 

The ICM mass profile is shown in figure [ij. We find 
Mgas = 2.56 ± 0.08 X 10" Mq within r = 1.7 Mpc or 
0.75 r2oo with a confidence level of 68%. 

We want to stress that we underestimate the central 
gas density in the centre due to neglecting the cooling flow 
excess. However, this underestimation at small radii does 
not play a significant role when we integrate the gas mass 
out to large radii. 

In fact, when we fit the surface brightness profile 
(figure |ll|) so that we only account for the inner 210 kpc, 
we find /3=0.49, rc=52kpc and Ueo — 5.55 x 10^^ cm^"^. 
By integrating equ. (^ with these parameters, we find a 
gas mass at r=0.1r2oo, which is 18% higher than the one 
using the global /3-model fit parameters. At r = 0.75r2oo, 
which is the radius up to which we can constrain the tem- 
perature and thus total mass profile we underestimate the 
gas mass by 1% due to neglecting the gas mass linked to 
the central gas excess with respect to our global /3-model 
, which is negligible. 
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Fig. 13. Abell 1835 total mass profile using the psf cor- 



rected temperature profile (see section 4.3). Error bars are 
at 68% confidence level(same display convention as figure 

HI). 



Fig. 14. Abell 1835 gas mass fraction profile. Error bars 
are at 68% confidence level (same display convention as 
figure HI). 



5.2. Cluster total mass 

Assuming spherical symmetry and hydrostatic equilib- 
rium, we calculate the gravitational mass of the cluster 
Abell 1835 with the hydrostatic equation : 



Mtot(< r) 



fcT 



dim 



dlnT 



d In r d\nr 



(7) 



where k is the Boltzmann constant, T the gas tempera- 
ture, G the gravitational constant, fj, the mean molecular 
weight of the gas (/i 0.6), rup the proton mass and rig 
the electronic density. Including equ. (^ in equ. (^, M can 
be written : 



Mtot(< r) 




(8) 



In order to account for the error bars in the temper- 
ature profile, we use a Monte-Carlo method described 
by Neumann & Bohringer (1995). This method allows to 
transform error bars of the temperature profile into error 
bars of the mass profile. This transformation is done by 
calculating randomly temperature profiles, which fit in the 
actual measured temperature profile and to subsequently 
calculate the mass profile corresponding to the randomly 
determined temperature profile. We determine 10000 ran- 
dom temperature profiles. From the corresponding 10000 
mass profiles we calculate the mean mass profile and the 
errors. The uncertainties are determined by taking into 
account the uncertainties of the temperature profile (90% 
confidence level) and the standard deviation of the mass 
at a given radius. The resulting errors are calculated so 
that they correspond to lu error bars (68% confidence 



level). The mass profile is shown in Fig.|l^. At a radius 
corresponding to 75% of the virial radius (or 1.7 Mpc) we 
obtain a total mass| of : Mtot = 1-23 ± 0.18 x 10^^ M© 
(68% confidence level). 

Assuming that Abell 1835 is isothermal with a temper- 
ature of 7.6 keV (see section iA), equ. (^) can be written 
as : 



Mtot(< r) 



3k(3 



T- 



(9) 



In this case taking into account error bars on T, /3 and rc 
parameters, the total mass of Abell 1835 is 1.00 ± 0.06 x 
10^^ Mq which is consistent with the mass obtained using 
the Monte-Carlo approach. This was expected since the 
temperature profile of Abell 1835 (see figure p^ ) seems to 
be isothermal at high radii. 

5.3. Gas mass fraction 

The gas mass fraction is simply the ratio of the ICM mass 
to the total mass. The profile of the gas mass fraction of 
Abell 1835 is shown in figure |l^. We do not display the gas 
mass fraction below 10% of the virial radius because of the 
cooling flow regio n an d our underestimate of the gas mass 
there (see section 4^). The errors bar are calculated with 
classical error propagation and are at a confidence level of 

1.5 



0.75r2oo, we find /gas/i^o = 0.207 ± 



68%. Within 

0.037 where h^o is the present time Hubble constant in 
units of 50km/s/Mpc. 



^ For a comparison, the total mass is 1.14 ± 0.14 x 10^^ Mq 
when we use the temperature profile without PSF correction 
(table §). 
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6. Discussion 

6.1. The gas mass fraction and implications for 
cosmology 

Abell 1835 is an ideal cluster to determine the gas mass 
profile with the hydrostatic approach, since the Xmm- 
newton image does not show indications of substructure 
at large scales]^, which might contaminate the mass deter- 
mination. 

Our determined gas mass fraction profile of Abell 1835, 
which we determine up to 0.75 r2oo is essentially flat. The 
mean value is about 21% and within 0.75 r2oo we obtain a 
gas mass fraction of 0.207 ± O.O37/I5Q . We can take this 
gas mass fraction as lower limit for the baryon fraction in 
the cluster because we do not take into account other con- 
tributions of baryons, such as stars or interstellar medium. 
Extrapolating this lower limit through out the universe 
(which was shown to be valid by White et al. (1993) 
and references therein), we can determine the cosmologi- 
cal matter density parameter Qm by combining our result 
with studies based on primordial nucleosynthesis, which 
give firm constraints on the density parameter of baryons 
in the universe. Recent results give ri^/igQ = 0.080 ±0.008 
(Buries et al. 2001). We can thus apply : 



gas 



Mu 



0.207 ±0.037/1 



-3/2 

50 



0.08 ±0.008/1 



50 



(10) 



If the ICM was the only source of baryons in clusters, above 
equation gives us ri,„ — 0.386 ± 0.108 (68% confidence 
level). However, since there exist other baryon sources in 
Abell 1835 such as galaxies, this implies that our result is 

— 1/2 

an upper limit with : < O.5/15Q . This result is in very 
good agreement with other recent results based on CMB 
measurements combined with redshift surveys (Efstathiou 
et al. 2001) and studies on distant supernovae (Perlmutter 
et al. 1999) which favor a low density universe with a 
cosmological constant. If we additionally knew the baryon 
fraction in galaxies up to high precision, we could put even 
tighter constraints on Qm- 

6.2. The temperature profile - comparison with other 
recent studies 

A lot of effort has been spent already for the determina- 
tion of temperature profiles in the ICM. The corresponding 
results vary from each study to another quite significantly. 
There have been reported declining temperature profiles, 
with a more or less universal shape based on AsCA and 
Beppo-SAX data (Markevitch et al. 1998, de Grandi & 
Molendi 2001). Other studies based on data from the same 
instruments, on the contrary, find flat temperature profiles 
or profiles with large dispersions (Irwin & Bregman 2000, 
White et al 2000). 



Schmidt et al. (2001) reported small scale substructure in 
the Chandra data in the cooling flow region, which does not 
influence our results on large scales. 



2 I 

0,0 



0.2 0.4 

Scaled radius (r/r_v) 



Fig. 15. Abell 1835 temperature profile and the compari- 
son with the universal profile derived by Allen et al. (2001) 
and extrapolated to 0.75 r2oo (dotted line) 



All these studies are limited to radii smaller than 
0.5r2oo due to instrumental limitations. Furthermore, 
studies based on AsCA and Beppo-SAX were subject to 
possible systematic uncertainties due to the large point 
spread function of the mirrors, which additionally can vary 
in size depending on energy. Chandra results presented 
by Allen, Schmidt & Fabian (2001) suggest that massive, 
relaxed clusters are isothermal up to r25oo and show a 
universal temperature decline in the center linked to the 
cooling flow. 

Xmm-newton with its better spatial resolution than 
AsCA and Beppo-SAX and larger effective area than com- 
parable observatories allows to a) better constrain the 
temperature profiles and b) to trace them out to larger 
radii. So far studies based on Xmm-NEWTON show either 
flat profiles (Arnaud et al. 2001a ; Tamura et al. 2001) 
for hot clusters, which is in agreement with the results by 
Allen et al. (2001) on Chandra data of clusters. A Xmm- 
newton study on the cool cluster Sersic 159-03, however, 
shows a dramatic temperature decline in the outer parts 
of the cluster (Kaastra et al. 2001). 

Our results suggest strongly a flat temperature profile 
with a decline towards the center linked to the cooling 
flow present in Abell 1835 (Peterson et al. 2001, Schmidt 
et al. 2001 and references therein). This profile fits well 
with the empirically found model by Allen et al. (2001), 
which is valid up to ^2500, which corresponds roughly to 
0.3 X r2oo and which is displayed in Fig |l^. However, an 
individual study undertaken by Schmidt et al. (2001) on 
Abell 1835 observed with Chandra shows a temperature 
profile which is increasing with radius. At radii larger than 
400 kpc Schmidt et al. (2001) find temperature estimates 
of about 12keV, which is inconsistent with our Xmm- 
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parts of the cluster, which agrees well with our derived 
temperature profile (see figure Hq) at large scales. 
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Fig. 16. Temperature profile of Abell 1835 corrected for 
the PSF (from table |[ crosses) comparing to the Chandra 
one (squares) obtained by Markevitch (2002). 



NEWTON results. The fact that we find good agreement 
between the individual spectral modeling of MOS2 and pn 
cameras gives us confidence in our results. 

While undertaking our study concerning the back- 
ground subtraction we observed that minor differences in 
the extraction of flares between background and source 
observation can lead to declining or increasing tempera- 
ture profiles. Doing monitoring on the background level we 
find that quiescent periods do not always have the same 
background level and have sometimes, even though rarely, 
intensities close or above the threshold criterion to extract 
flares. In order to screen the data properly for background 
variations and to determine the level of the particle (in- 
strumental) background accurately it is important to look 
at parts of the observation where no astrophysical source 
emission can be detected like the high energy part of the 
spectrum, where only particle background can be seen 
(the particle background has a hard spectrum and the 
effective area of Xmm-newton at energies above lOkeV 
for X-ray photons is relatively small). It is important to 
compare and to apply exactly the same flare rejection for 
both the source and the background observation to avoid 
wrong background subtraction, which can bias the tem- 
perature estimates. Observing a temperature profile which 
increases with radius can be an indication of a slightly 
higher particle background in the source observation than 
in the background observation. 

The presence of remaining solar flares in the anal- 
ysis of Schmidt et al. (2001) was recently discussed by 
Markevitch (2002), who reanalyzed the Chandra data of 
Abell 1835. Markevitch (2002) found that a temperature 
profile analysis based on a more thorough flare treatment 
leads to a basically flat temperature profile in the outer 



6.3. Global temperature estimates 

There is some dispersion observable between past tem- 
perature determinations of Abell 1835 based on AsCA- 
data. Allen & Fabian (1998) found a temperature of 
9.8j;^ 3keV, based on AsCA-SiS and GiS-data, including 
a multiphase cooling flow model, which is now known 
not to be an accurate description for the cooling flow of 
Abell 1835 (Peterson et al. 2001). White (2000) only used 
AsCA-GiS data and obtained an overall temperature of 
kT=8.5lQ5 keV when applying a multiphase cooling flow 
model and kT=8.2 ± 0.5 keV when using a single temper- 
ature model. Both results are consistent within the error 
bars with our results of kT=7.6±0.4keV. Furthermore, 
Mushotzky & Scharf (1997) found a temperature estimate 
of kT=8.2 ± 0.5 keV, which is again comparable with our 
results and the ones of White (2000). The results of Allen 
& Fabian (1998) agree within the error bars with other 
results based on AsCA, but are higher than our mean 
temperature. While our temperature results are somewhat 
lower than the results of AsCA, the Chandra results for 
the temperature of Abell 1835 by Schmidt et al. (2001) 
are above the results based on AsCA. 

7. Conclusion 

Xmm-newton data allow us to measure the temperature 
profile of Abell 1835 up to 0.75 r2oo- In order to determine 
the temperature profile of the cluster accurately we apply 
a method, which corrects for the various kinds of back- 
ground and a photon weighting method which allows to 
correct for vignetting effects. We correct the temperature 
profile for PSF-effects and see that the resulting profile is 
not affected by the PSF at large radii. 

Our resulting PSF corrected temperature profile is con- 
sistent with being flat in the outer parts. The mean cluster 
temperature is 7.6 keV outside the cooling flow region. In 
the central parts, below a radius of l'(300kpc) we see a 
temperature decline linked to the cooling flow of the clus- 
ter. We fit a /3-model to the outer radii of the surface 
brightness profile of Abell 1835 and obtain the best fit re- 
sults of /3 = 0.704 ± 0.005 and r^ = 202.3 ± 7.1 kpc. We 
use these best fit parameters as input for the hydrostatic 
approach to calculate the mass profile of the cluster and 
find a mass at 0.75 raoo of 1.23±0.18x lO^^Mg. The corre- 
sponding gas mass fraction, which is constant with radius 
at radii r > 0.1r2oo, is with 20.7%±3.7 comparable with 
the results found in other studies. 

Our global temperature estimate for Abell 1835 based 
on Xmm-newton is somewhat low when compared to 
other results, but agree in the error bars with other studies 
based on AsCA-data. We assume that these discrepencies 
are linked to small remaining uncertainties in the calibra- 
tion of Xmm-newton , and on which is still being worked 
on. 
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